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[57] ABSTRACT 

Apparatus and methods are disclosed for the detection and 
imaging of ultrasonic harmonic contrast agents. The har- 
monic echo effect is detected through alternate polarity 
acquisition of harmonic contrast agent effects, which pro- 
vides the benefits of suppressing the harmonic components 
of the transmitted signal while eliminating clutter. 

4 Claims, 5 Drawing Sheets 
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ULTRASONIC DIAGNOSTIC IMAGING 
WITH HARMONIC CONTRAST AGENTS 

This invention relates to ultrasonic diagnosis and imag- 
ing of the body with ultrasonic contrast agents and, in 
particular, to new methods and apparatus for ultrasonically 
detecting and imaging with contrast agents. 

Ultrasonic diagnostic imaging systems are capable of 
imaging and measuring the physiology within the body in a 
completely noninvasive manner. Ultrasonic waves are trans- 
mitted into the body from the surface of the skin and are 
reflected from tissue and cells within the body. The reflected 
echoes are received by an ultrasonic transducer and pro- 
cessed to produce an image or measurement of blood flow. 
Diagnosis is thereby possible with no intervention into the 
body of the patient 

However materials known as ultrasonic contrast agents 
can be introduced into the body to enhance ultrasonic 
diagnosis. Contrast agents are substances which will 
strongly interact with ultrasonic waves, returning echoes 
which may be clearly distinguished from those returned by 
blood and tissue. One class of substances which has been 
found to be especially useful as an ultrasonic contrast agent 
is gases, in the form of tiny bubbles called microbubbles. 
Microbubbles present a significant acoustic impedance mis- 
match in comparison to tissue and fluids, and nonlinear 
behavior in certain acoustic fields which is readily detectable 
through special ultrasonic processing. In order to infuse 
bubbles into the body so that they will survive passage 
through the pulmonary system and circulate throughout the 
vascular system, gases have been stabilized in solutions in 
the form of tiny microbubbles. Microbubble contrast agents 
are useful for imaging the body's vascular system, for 
instance, as the contrast agent can be injected into the 
bloodstream and will pass through the veins and arteries of 
the body with the blood supply until filtered from the blood 
stream in the lungs, kidneys and liver. 

One property of microbubble contrast agents currently 
under investigation is harmonic response. These harmonic 
contrast agents exhibit significant, detectable responses at 
frequencies which are harmonics of the transmitted ultra- 
sonic frequency. This property is useful for clutter rejection 
of the received signals. When the transmitted frequency 
band is used as the received frequency band, echoes will be 
returned from the microbubbles, but also from surrounding 
tissue, the latter comprising clutter in the received echo 
signals. But with harmonic contrast agents, reception occurs 
at harmonic frequencies, where fundamental band clutter 
from tissue is ignored. Since tissue generally reflects very 
minimal harmonic components, the received harmonic band 
enables the microbubble echoes to be received with a high 
signal to noise ratio. 

In accordance with the principles of present invention, a 
technique is provided for the detection and imaging of 
harmonic ultrasonic contrast agents. The harmonic contrast 
agent is insonified by alternate polarity transmitted pulses, 
and the echo signals received from the transmitted pulses are 
combined. The result is a suppression of harmonic compo- 
nents of the transmitted ultrasonic waves and the elimination 
of clutter. 

In the drawings: 

FIG. 1 illustrates in block diagram form apparatus con- 
structed in accordance with the present invention which 
provides performance advantages for harmonic contrast 
agent detection; 

FIGS. 2 and 3 illustrate passband characteristics used to 
explain the performance of the embodiment of FIG. 1; 
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FIGS. 4a and 4b illustrates the alternate polarity pulsing 
of harmonic contrast agents; and 

FIGS. 5a-Sc illustrate nonlinear response waveforms 
produced by alternate polarity acquisition of contrast agent 
5 echoes. 

Referring to FIG. 1, an ultrasonic diagnostic system for 
use with harmonic contrast agents in accordance with the 
present invention is shown in block diagram form. In this 
system an array transducer 112 of a probe 110 transmits 
ultrasonic energy and receives echoes returned in response 
10 to this transmission. The response characteristic of the 
transducer can exhibit two passbands, one around the central 
transmit frequency and another about the center of the 
received passband. For imaging harmonic contrast agents, a 
broadband transducer having a passband encompassing both 
15 the transmit and receive passbands is preferred. The" trans- 
ducer may be manufactured and tuned to exhibit a response 
characteristic as shown in FIG. 2, in which the lower hump 
60 of the response characteristic is centered about the center 
transmit frequency f„ and the upper hump 62 is centered 
20 about the center frequency f r of the response passband. The 
transducer response characteristic of FIG. 3 is preferred, 
however, as the single dominant characteristic 64 allows the 
probe to be suitable for both harmonic contrast imaging and 
imaging without harmonic contrast agents. The characteris- 
25 tic 64 encompasses the central transmit frequency f„ and 
also the harmonic receive passband bounded between fre- 
quencies f L and f c , and centered about frequency f r A typical 
harmonic contrast agent can have a response such that 
transmission about a central transmit frequency of 1.7 MHz 
30 will result in harmonic returning echo signals about a 
frequency of 3.4 MHz. A bandwidth characteristic 64 of 
approximately 2 MHz would be suitable for these harmonic 
frequencies. 

In FIG. 1 a central controller 120 provides a control 
35 signal i„ to a transmit frequency control circuit or pulser 121 
to control the center frequency and time of transmission of 
the transmitted ultrasonic energy. The transmit frequency 
control circuit pulses the elements of the transducer array 
112 by means of a transmit/receive switch 114. 
40 Echoes received by the transducer array 112 are coupled 
through the T/R switch 114 and digitized by analog to digital 
converters 115. The sampling frequency f, of the A/D 
converters 115 is controlled by the central controller. The 
desired sampling rate dictated by sampling theory is at least 
45 twice the highest frequency f c of the received passband and, 
for the preceding exemplary frequencies, might be on the 
order of at least 8 MHz. Sampling rates higher than the 
minimum requirement are also desirable. 

The echo signal samples from the individual transducer 
50 elements are delayed and summed by a beamformer 116 to 
form coherent echo signals. The digital coherent echo sig- 
nals are then filtered by a digital filter 118. In this 
embodiment, the transmit frequency f fr is not tied to the 
receiver, and hence the receiver is free to receive a band of 
55 frequencies which is separate from the transmitted band. The 
digital filter 118 bandpass filters the signals in the passband 
bounded by frequencies and f 0 in FIG. 3, and can also shift 
the frequency band to a lower or baseband frequency range. 
The digital filter could be a filter with a 1 MHz passband and 
60 a center frequency of 3.4 MHz in the above example. A 
preferred digital filter is a parallel arrangement of serially 
coupled multipliers and accumulators. This arrangement is 
controlled by the central controller 120, which provides 
multiplier weights and decimation control which control the 
65 characteristics of the digital filter. Preferably the arrange- 
ment is controlled to operate as a finite impulse response 
(FIR) filter, and performs both filtering and decimation. 
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Filtered echo signals from tissue, generally filtered by a 
passband centered about or demodulated from the transmit 
frequency, are coupled to a B mode processor 37 for 
conventional B mode processing. Filtered echo signals of 
the harmonic contrast agent passband are coupled to a 5 
contrast signal detector 128 which performs pulse to pulse 
summation or integration of temporally discrete echoes from 
a given spatial location, amplitude or envelope detects the 
combined signals. Simple two pulse summation of the form 
P!+P 2 may be employed where P 1 represents the echoes 10 
received following one pulse and P 2 represents the echoes 
received following another pulse. The combination of ech- 
oes from consecutive pulses may, if desired, be performed 
before the digital filter 118 rather than after, the decision 
being a matter of choice of system design. 15 

The filtered echo signals from the digital filter 118 are 
also coupled to a Doppler processor 130 for conventional 
Doppler processing to produce velocity and power Doppler 
signals. The outputs of these processors are coupled to a 3D 
image rendering processor 132 for the rendering of three 20 
dimensional images, which are stored in a 3D image 
memory 134. Three dimensional rendering may be per- 
formed as described in U.S. patent application Ser. No. 
08/638,710, and in U.S. Pat. Nos. 5,474,073 and 5,485,842, 
the latter two patents illustrating three dimensional power 25 
Doppler ultrasonic imaging techniques. The signals from the 
contrast signal detector 128, the processors 37 and 130, and 
the three dimensional image signals are coupled to a video 
processor 140 where they may be selected for display on an 
image display 50 as dictated by user selection. The video 30 
processor preferably includes persistence processing, 
whereby momentary intensity peaks of detected contrast 
agents can be sustained in the image. One technique for 
providing persistence is through frame averaging, whereby 
new image frames are combined with previous frame infor- 35 
mation on a spatial basis. The combination can be done by 
weighting the contributions of the old and new frame 
information and the frame information can be combined in 
a recursive manner; that is, old frame information is fed back 
for combining with new frame information. A preferred 40 
persistence technique is the fast attack, slow decay technique 
described in U.S. Pat. No. 5,215,094, which can be applied 
to both Doppler and contrast agent images. 

The apparatus of FIG. 1 performs alternate polarity pulse 
transmission as illustrated in FIGS. 4a and 4b. In the first 45 
transmission of FIG. 4a, the central controller 120 provides 
a first polarity control signal f^ to the pulser 121, which 
drives the transducer elements 112 to transmit a first polarity 
pulse 20. For the second transmission of FIG. 4b, the central 
controller 120 provides a second polarity control signal f^, 50 
to the pulser 121, which drives the transducer elements 112 
to transmit a second polarity pulse 22. 

The echoes received from microbubbles in response to 
these alternate polarity transmissions are shown in FIGS. 5a 
and 5*. FIG. 5a illustrates an echo waveform 310 received 55 
from the first pulsing of a microbubble contrast agent. The 
nonuniform amplitudes on either side of the zero reference 
level illustrate nonlinear reflexive action of microbubbles in 
the presence of acoustic waves, as the microbubbles non- 
linearly compress and expand. The echo waveform of 310 60 
FIG. 5a results from transmission of an ultrasonic pulse 
exhibiting a first polarity. 

Following transmission of the ultrasonic pulse exhibiting 
the opposite polarity, the echo waveform 312 of FIG. 5b 
results. This waveform is similarly nonlinear, but out of 65 
phase with the first waveform due to the change in pulse 
polarity. When the two waveforms are combined, a har- 



monic response is obtained, as shown in FIG. 5c. The highly 
nonlinear waveform of FIG. 5c is readily detected, causing 
the system to become highly sensitive to the contrast agent 
which produced the nonlinear echo responses. 

A mathematical analysis of this effect and response is as 
follows. To detect the harmonic response of microbubbles, 
the harmonic component in the incident pressure wave must 
be suppressed. Based on the analytical solution of the 
dynamic motion of microbubbles, the primary component of 
the backscattering pressure magnitude is linearly propor- 
tional to the incident pressure and the harmonic component 
is quadratically proportional to the incident pressure p, or 
pXceOccp,- and p^cojapj 2 . Thus, neglecting the higher order 
terms, one may write the backscattering pressure magnitude 
p B ((Q) from a microbubble in a generic form 



(1) 



where k t and are parametrically related to the acoustic 
properties of the microbubble such as size, viscosity, surface 
tension, ambient pressure, etc. 

Now assume that the microbubble is excited by two 
narrow band signals at different times but with the same 
magnitude p and at the same frequency co, but with opposite 
polarity: p ; i=p cos cot and p,- 2 =-p cos cot Then the back- 
scattered pressure wave from Pn=p cos cot is 

P*i(a>>')=*i(®,f)/H-*2(a>,')p 2 (2) 

and from p i2 =-p cos cot is 

p J2 (co,r+5r)=t 1 ((»,(+5() P +jt 2 (a),t+50p 2 (3) 

Then the total backscattered pressure magnitude may be 
obtained by summing Equations (2) and (3), 



S = P Bl + pjB = (tl(0),t)-tl(0),f+8«))p + (fe(CO,f) + fe(<B,f + 8f))F ! W 

Equation (4) shows that the primary component is elimi- 
nated if k^co) and lc^co) do not change substantially in the 
time duration 6t, where 8t is small. 

Assume the backscattering from microbubbles is quasi- 
stationary over T, where T is the pulse repetition interval. 
Therefore, the average nonlinear acoustic properties are not 
changed over time T, or 

£{fc 1 (co,»)}=B{i 1 (to,H-7)} 



E{i 2 (G»)}^{fc 2 (t 0 ,(+D}. 

The relationship of Equation (4) will hold by summing the 
pulse echoes from two pulses which are time-diverse in T. 
The quasi-stationary assumption is valid for slow perfused 
flow, such as myocardial perfusion. 

When the bandwidth of the incident pressure wave is 
wide, the wideband excitation wave P(t) may be represented 
by a Fourier series 

P(t) = X ACcUj^osov 



Thus the backscattered pressure magnitude of the 
microbubbles from P(t) may be written as 
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and the backscattered pressure magnitude of the 
microbubbles from -P(t) may be written as 



Accumulating the partial s 
results in 



m of consecutive pairs of echoes 



Pbi = - Z fci(a>;yt(05;) + 1 feCmOA^wi) 
Summing Equations (5) and (6), one may obtain 

S=pBi+PBi = 2-Lk 2 (,<s> i yf(a i ) 0) 

Again, the harmonic component is extracted and the primary 
component is eliminated. 

Let us assume the nonlinearity in tissue is negligible. 
Since the backscattered pressure in a linear medium is 
linearly proportional to the incident pressure wave, the 
polarity of the backscattered wave will be changed as the 
polarity of the incident pressure wave is changed. Assuming 
the tissue is relatively stationary during the period of two 
consecutive pulses, summing the pulse echoes from con- 
secutive pulses with opposite polarity will cancel the echo 
response from tissue. Thus, tissue clutter will be suppressed. 

The concept of summing the pulse echoes from two 
pulses of opposite polarity may be generalized into process- 
ing echoes from multiple pulses with alternate polarity to 
maximize the sensitivity and ininimize the variance, assum- 
ing the tissue is stationary during the pulsing interval. Let 
the pulse sequence be 



P={p -p p- p p- p "-1 
and the pulse echoes be 

E={E 1 E 2 E 3 E 4 E 5 E 6 ' 



p} 



S= £ £ f + E /+ i = 2(n-l)ri 2 (0)i)A 2 (a)i) 
j=l i 

* What is claimed is: 

1. A method of ultrasonically detecting the ultrasonic 
10 response of an ultrasonic contrast agent comprising the steps 

of: 

transmitting a first ultrasonic pulse to said ultrasonic 
contrast agent to cause a first harmonic response; 
15 transmitting a second ultrasonic pulse of a different polar- 
ity than said first ultrasonic pulse to said harmonic 
contrast agent to cause a second harmonic response; 
detecting said first and second harmonic responses; and 
20 combining said first and second harmonic responses. 

2. The method of claim 1, wherein said step of combining 
comprises summing said first and second harmonic 
responses. 

3. The method of claim 1, wherein said step of combining 
25 comprises integrating said first and second harmonic 

responses. 

4. The method of claim 1, wherein said transmitting step 
comprises transmitting pulses which exhibit a pulse energy 

30 which is within a range which causes microbubbles of said 
ultrasonic contrast agent to oscillate without substantial 
microbubble destruction. 



